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. Abstract. We studied the formation process of star clusters using high-resolution A^-body/smoothed 

Oh- particle hydrodynamcs simulations of colliding galaxies. The total number of particles is 1.2 x 10 

q , for our high resolution run. The gravitational softening is 5 pc and we allow gas to cool down to 

!— i ~ 10 K. During the first encounter of the collision, a giant filament consists of cold and dense 

' gas found between the progenitors by shock compression. A vigorous starburst took place in 

. the filament, resulting in the formation of star clusters. The mass of these star clusters ranges 
from 10 5 " 8 Mq. These star clusters formed hierarchically: at first small star clusters formed, 
and then they merged via gravity, resulting in larger star clusters. 
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1. Introduction 



Merging galaxies contain many young star clusters (e.g., Whitmore 2003). These star 
clusters could potentially evolve into the present day metal-rich globular clusters and so 
they are widely accepted to be a good candidate for globular cluster progenitors. Galaxy- 
galaxy merger is therefore considered as one of the formation channels of globular clusters 
/S ; (Schweizer 1987). 

There have been a number of numerical studies of merging galaxies. There have been, 
" ' however, only a few numerical studies of star cluster formation in merging galaxies even 
though it has been shown that resolving a cloudy /multiphase interstellar medium (ISM) 
and/or clustered star formation can have important consequences for the formation his- 
tory of early-type galaxies (e.g. Bois et al. 2010). Some of the existing studies adopted 
sub-grid models of star cluster formation (e.g., Bekki & Couch 2001; Li ct al. 2004), while 
more recently their formation has been captured directly (e.g. using the sticky particle 
method in Baurnaud et al. 2008). Here we report the result of merger simulations that 
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capture the multiphase nature of the ISM and include realistic models of star formation 
and feedback. 



2. Method 

We prepared two identical progenitor galaxies and then let them merge from a parabolic 
and coplanar configuration. The mass of components in one progenitor galaxy is 10 11 Mq 
for the dark matter halo, 4.7 x 10 9 M for the stellar disk and 1.8 x 10 9 Mq for the gas 
disk. The galaxies are modeled using both TV-body and smoothed particle hydrodynamics 
(SPH) particles. We employed 1.2 x 10 s particles for the two progenitor galaxies for the 
finest runs where the corresponding mass of each particle is 1.9 x 10 3 Mq for both TV-body 
and SPH particles. 

Numerical simulations were performed by our TV-body/SPH code, ASURA. We adopted 
the FAST scheme (Saitoh & Makino 2010) that accelerates the time integration of a self- 
gravitating fluid by asynchronously integrating gravity and hydrodynamical interactions. 
In addition, we used a time-step limiter for the time-integration of SPH particles with 
individual time-steps which enforces the differences of time steps in neighboring SPH 
particles to be small enough so that the SPH particles can correctly evolve under strong 
shocks (Saitoh & Makino 2009). 

The high mass resolution allows us to adopt the ISM model with the wide temperature 
range down to 10 K and realistic conditions of star formation (p > 100 cm -3 , T < 100 K 
and V -v < 0, where p, T and v indicate density, temperature and velocity, respectively). 
The feedback from type II supernovae (SNe) was also taken into account. These models 
are the same as those in Saitoh et al. (2008; 2009). 

3. Results 

At the first encounter (after ~ 420 Myr from the beginning of the simulations) , strong 
shocks took place. These shocks induced the formation of a large filament of cold and 
dense gas between the two galaxies. This filament quickly cooled and fragmented to 
a number of small, high-density clumps, and a large-scale starburst took place in the 
filament (see Saitoh et al. 2009). This starburst did not take place in previous studies in 
which gas cooling was inhibited below 10 4 K. This starburst continued for several 10 Myr 
and was then quenched by the SN feedback. A number of star clusters formed during 
this starburst. 

Figure [1] shows the snapshots of the early phase of the star cluster formation where 
we can see that a number of small star clusters formed along the gas filament. The late 
phase of the star cluster formation (i ^ 440 Myr) was mainly driven by mergers of star 
clusters without gas, since gas in the filament was blown out by SNe. Hence star clusters 
grew hierarchically through gravitational mergers of smaller star clusters. There is no 
star cluster with the mass ^ 10 7 Mq which forms from an instantaneous collapse of a 
single large cloud. At the final phase (after 40 Myr from the start of the starburst), we 
found several tens of star clusters between two progenitor galaxies (see figure [2]) . This is 
because they were formed only from the gas in the shock-induced filament. 

Figure [3] shows cumulative mass functions of star cluster systems. The mass of star 
clusters ranges from 10 5-8 Mq. The mass function for masses larger than 10 7 Mq is 
almost independent of mass and spatial resolutions. The slope of this high mass region is 
close to a power law function with an index of -2, which is consistent with observations 
(e.g., Whitmorc et al. 1999). Since the formation of star clusters is driven by mergers, 
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Figure 1. Distributions of gas and stars in the early phase of the star cluster formation 
(t = 426 Myr). Main panel shows gas density map of the mid-plane in 12 kpc x 12 kpc re- 
gions, while the bottom-left inset displays four times larger scale than that of the main panel. 
Closeup views of star-cluster forming regions show surface stellar density maps with contours 
of surface gas density in 500 pc x 500 pc and 100 pc x 100 pc. 

the shape of the mass function becomes a power law reflecting the scale free nature of 
gravity (Elmegreen 2006). 

4. Summary 

We performed high resolution simulations of merging galaxies that capture the mul- 
tiphase ISM and the realistic modelings of star formation and of SN feedback. In the 
simulations, we found that a number of star clusters formed during the starburst at the 
first encounter. These star clusters grew via hierarchical mergers. 
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Figure 3. Cumulative mass functions of star cluster systems at t = 460 Myr. Thick and thin 
curves indicate cumulative mass functions of star clusters in runs with e = 5 pc and 20 pc, 
respectively. Solid, dashed, dotted curves represent high, middle, low mass resolution runs, 
respectively. 
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